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We consider the Casimir interaction, mediated by massless fermions, between a spherical defect 
and a flat potential barrier, assuming hard (bag-type) boundary conditions at both the barrier 
and the surface of the sphere. The computation of the quantum interaction energy is carried out 
using the multiple scattering approach, adapted here to the setup in question. We find an exact 
integral formula for the energy, from which we extract both the large and short distance asymptotic 
behaviour. At large distance the fermionic contribution is found to scale as L~^, in contrast to 
that of electromagnetic vacuum fluctuations that, assuming perfectly conducting boundaries, scales 
as L~^. At short distance, we compute the leading and sub-leading contribution to the vacuum 
energy. The leading one coincides with what it is expected from the proximity force approximation, 
while the sub-leading term gives, contrary to the electromagnetic case, a positive correction to the 
proximity force result. 
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Introduction. In 1948 Casimir predicted that two flat, 
and perfectly conducting plates would modify the vac¬ 
uum fluctuations of the electromagnetic field and in¬ 
duce an observable force between them [I|. This set-up 
was difficult to unambiguously test experimentally at the 
time 0 , and the first conclusive observation came many 
years later, in fact, for the different configuration of a 
sphere near a plate Q. It took some more time to con¬ 
clusively observe the Casimir force between two plates, 
which was, eventually, achieved a few years later (See 
Ref. Q for review). 

Much of the work that followed Casimir’s conclusion 
focused on the electromagnetic vacuum fluctuations and 
different geometries (See Refs. @). However, the par¬ 
allel plates setup was just one specific example of the 
more generic deformations in the quantum vacuum that 
boundaries can produce, irrespectively of the nature of 
the quantum fields. The key feature is the presence of 
massless (or quasi-massless) quanta that induce long- 
ranged correlations, suggesting that quantum vacuum ef¬ 
fects may be relevant in fermionic environments. The 
analogous of the electromagnetic Casimir effect should, 
in fact, occur in condensed matter systems, like in quan¬ 
tum liquids @ , or when long-range correlations exist due 
to Goldstone modes of a broken continuous symmetry, as 
in superfiuids 0, [1| . 

One particularly exciting example is associated with 
the presence of defects in quantum fermi liquids, and with 
the possibility that high precision emerimental manipu¬ 
lation of ultra-cold atomic systems , where defects can 
be controlled in a variety of ways, may offer novel tests of 
quantum vacuum energy effects and, for instance, provide 
new constraints on hypothetical sub-micron interactions. 


Other physical setups relevant to the fermion Casimir 
effect include, for instance, carbon nanotubes Elm, 
graphene MM , nuclear structures M , neutron star 
crusts [1^ (See Refs, [l^ for a longer fist of examples 
and references). 

The simplest setups analysed have, so far, ignored the 
structure of the defects and focused on one-dimensional 
boson and fermi systems with defects treated as delta- 
functions in the adiabatic approximation (See, for exam¬ 
ple, Refs. [iS[il,[i3)- In Refs. [ 1 ^, it was found that, for 
a generic (interacting or noninteracting) fermionic back¬ 
ground, the Casimir force between the impurities oscil¬ 
lates as a function of the separation. The similar problem 
of calculating the Casimir interaction between two scat- 
terers immersed in a one-dimensional massless fermionic 
background has also been analysed in Ref. M , using a 
force operator approach, with the defects modelled by 
two delta-function potentials. A similar oscillatory be¬ 
haviour has been found too, however, as a function of 
relative polarization of the two scatterers, while the de¬ 
pendence on the distance turned out to be monotonic 
(for fixed polarizations). 

A fuller understanding of the Casimir effect mediated 
by fermions should include the structure of the defect and 
extend to higher dimensionality. This class of problems 
has been analysed, for instance, in Refs, [l^, where the 
properties of systems containing one or more fermionic 
bubbles (almost spherical defects immersed in a homoge¬ 
nous fermionic environment) have been discussed. The 
Casimir energy of a system composed of two, three, and 
four spheres in a fermionic background has been analysed 
in Ref. [l8l| in the semiclassical approximation and led to 
the suggestion that in many-body Casimir interactions. 
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the two-body ones dominate at small separation. 

In the present work, we wish to consider the related 
problem of computing the Casimir interaction between 
a defect and a potential barrier, encapsulating a mass¬ 
less fermionic system. Following Ref. [^, we model the 
defect as a spherical bubble, and both the bubble and 
the barrier as hard walls. The geometry of the system is 
that of a sphere of radius R, centered at the origin, close 
to a plate of surface area H x H located at z = L, at 
a distance d from the sphere. Since we are considering 
the case of a defect whose size is much smaller than the 
surface area of the wall, we shall assume that R H. 

Despite the simplicity of the setup, the problem of 
studying the Casimir interaction between a sphere and a 
wall is not straightforward and has been subject of many 
analyses that have focused on the scalar and electromag¬ 
netic quantum vacuum fluctuations [20I - [T7[ [^ . While 
the leading term in this interaction at short distance can 
be obtained easily using the proximity force approxima¬ 
tion [l^, going beyond is difficult. 

The increasing precision at which Casimir force mea¬ 
surements can be performed called for more accurate 
computations and for efficient ways to go beyond the 
proximity force approximation. An especially advanta¬ 
geous one is based on the multiple scattering approach 
and has offered a systematic way to compute the Casimir 
energy between two compact objects of arbitrary shape 
[20I - [T7I . [ 2 ^. Results have covered a variety of cases, 
mainly for the scalar and electromagnetic fields (See, 
for instance, Refs. [23l - l^ or Ref. [ 2 ^ for review). The 
beauty of the multiple scattering approach is that it is 
physically transparent, since the interaction energy is 
expressed in terms of a multiple scattering expansion 
(waves that scatter back and forth between the two ob¬ 
jects that are interacting). It is also of straightforward 
numerical implementation, and it allows for a systematic 
way to extract sub-leading corrections to proximity force 
results. 

For the sphere-plate configurations, the multiple scat¬ 
tering approach has been adopted to obtain the scalar 
Casimir energy beyond the proximity force approxima¬ 
tion has been worked out in Refs. [^,|^. Results can be 
summarised in the formula Ri 1 -|- a^^d/R, 

where X and Y represent the boundary conditions im¬ 
posed at the plate (X = Dirichlet (D), Neumann (N), 
Robin (R)) and at the sphere (Y = D, N, R) , respec¬ 
tively, and = p^'^R/d?. The numerical coeffi¬ 

cients are pDD/NR = -Tr^/lddO, p^D/DR ^ _77r3/ii520, 
^DD/ND ^ ^dr ^ gQ( 3 ^ _ 2)/(77r2) and 

^NR _ 1/3 _|_ 20(3a — 2)/7r^ (a is the Robin parame¬ 
ter). For the electromagnetic field, the Casimir energy 
beyond the proximity force approximation has been com¬ 
puted in Refs. for perfectly conducting boundary 

conditions, leading to k, R/{J2Qd?')(\ -|- (1/3 — 

20/7r2)d/i?). 

Fermion TGTG formula. While the multiple scatter¬ 
ing approach has been systematised for the electromag¬ 
netic and scalar Casimir effect, only limited attention has 


been paid to the fermionic case. The possibility of high 
precision Casimir effect experiments in fermionic envi¬ 
ronments (for instance, in ultra-cold atomic systems Q), 
provides a natural motivation to carry out more precise 
computations, beyond proximity force results. 

As shown in Refs. [22 - 1^ . the multiple scattering ap¬ 
proach allows to express the Casimir energy in terms of 
transition matrices (associated to the interacting bod¬ 
ies) and the propagators, thus offering a precise com¬ 
putational prescription, and, in the following, we adapt 
this approach to our case. The formal expression for the 
fermion Casimir energy takes the usual form 

i^Cas = -^y d^lndet (I - N), (1) 


where N = and ^ is the imaginary frequency. 

The matrices T®, i = 1, 2 are the transition matrices 
associated to the boundaries (the sphere and the wall 
in our case) and G®-^ represent the translation matrices. 
Formally, aside for the change in the overall sign, in the 
expression above (also called TGTG-formula after the 
work [13), the nature of the quantum fields is encoded 
in the matrices inside the determinant. The matrices 
T® are related to the scattering matrix of object i and 
can be computed by matching the boundary conditions 
imposed on object i. The translation matrices G®-^ relate 
the basis of wave functions of object i to the basis of 
wave functions of object j. Expression o is valid at a 
formal level and deriving an explicit expression reduces to 
calculating the above matrices and taking determinant. 
Here we will follow the procedure outlined in Ref. [13, 
where a prescription to compute the matrices T® and G®-^ 
using the mode-summation approach has been developed. 

The fermions in our computation are massless spin 
1/2-fermionic fields "0, satisfying the Dirac equation, 
= 0, where = 9^+F^, and F^ is the spin con¬ 
nection. On the boundaries of the sphere {B = sphere) 


and the wall {B = wall), we impose (1 -I- 


= 0 , 


with being the unit outward normal vector. 

In order to match the boundary conditions and find the 
transition matrices, one needs to express the solutions in 
terms of a spherical and plane wave basis, respectively 
for the sphere and the plate. Explicit forms are known 
for spherical waves (See Ref. HI): 




(±),* 

jm,l 




(±).* 

jm,2 






( 2 ) 


where each mode is characterized by the quantum num¬ 
bers (j, m), with j = 1/2, 3/2, 5/2, ... and m = —j, —j + 
1,..., j — 1, j. The superscripts (+) and (—) indicate, re¬ 
spectively, the positive and negative energy modes, Dum, 
represent the spherical harmonic spinors (see Ref. [33j), 
and k = w/c. For convenience, we have adopted the 
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notation * = reg (* = out) for regular (outgoing) waves, 




with , C°'^* = TTp'^i/2. The plane waves 

can be parametrized in terms of the momenta perpendic¬ 
ular to the plate, kj^ = (fci, /C 2 ), 

_ j^^),*^ikix+ik2y-isgn,y/k'^-k\z^iLjt 

where kj_ = ^/kl + fc|, sgn^eg = 1, sgu^^^t = -1, and 

I 1 \ 


.(±).* _ 
"^kx.l “ 


A 


(±).* _ 

k, ,2 “ 


Tsgn, 

± 


± 


1 

0 _ 

y/k'^ — k"^ 

ki + 1^2 

k 

0 

1 

fci — ik2 


( 5 ) 


y±sgn* 


y/k'^ — k‘j 


The matching procedure is tedious but straightforward. 
It consists in expressing one set of waves in terms of the 
other, with the coefficients of this transformation defin¬ 
ing the transition matrices. Imposing the boundary con¬ 
ditions, and solving the resulting equation for the tran¬ 
sitions coefficients allows one to find and (Details 
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for this and all the other calculations are reported in ap¬ 
pendix). These are block-diagonal matrices in (j, m) and 
kj^ respectively. The (j, m)-block of is a diagonal 2x2 
matrix of the form 




T, 


(±) 


0 


0 


T, 


(±) 


where 


rpi^) _ 
j-m 


Ij (kR) =F Hj+i (kR) 

Kj{KR)±iKj+i{KR)' 

with K = ik. The kj^-block of is 




cos 9k sinSfeC 
— sin0fce®‘^'= cos 6 k 


( 6 ) 


( 7 ) 


where dk and (pk are defined so that y/k'^ — k\ = k cos 9k , 
fci = fc_L cos Pk and fc 2 = fc_L sin pk ■ 

The translation matrices and G^^ are dehned by 
the relations 


''C!r 


V'SrVx.^)' 





. 12 ,(±) 
F • 1 

jm,kx 




^(±).reg 


(x,a;) 


=H‘‘ 


(Pk 


.21.(±) 


VS!r(x',u;)^ 


( 8 ) 

where x' = x — L, L = (0,0, L). The computation can 
be performed following the idea introduced in [s^ for 
the computation of the corresponding translation matri¬ 
ces for scalar and electromagnetic fields. We find the 
following result : 




(j -b m)P'^_P {cos9k) {cos9k)e 

I i{j - m + 1)^™ (cosdfe) {cos6k)e 

J' 2 •?"*” 2 


—i^k 


(9) 


G; 


2i.(±) I U-m)\ jU + m)Pj^_,^cos 9 k) pj - m + l)P^^,^cos 9 k) 

" 2iJ2Y (j + m)!fc^fc 2 _fc^® 


P™t"(cos6»fe)e*‘^'= 
2 


-iiP”(t"(cos6»fe)e*^'' 
2 


( 10 ) 


Combining the above expressions (using some identities for the associated Legendre functions P;'"(z), noticing that 
the argument in the determinant is real, and the contribution from the positive energy modes and negative energy 
modes are the same), one can obtain the Casimir energy written as 


PCas 


d^lndet (I-M) 


( 11 ) 
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where 


M- , 


w / (j - m)!(j' - m)\ f _0_\ 

2 Y (j + m)!(j' + m)! q -Tj+7 

(j + m)P'^_ 1 ^ (cosh 0) ^ (cosh 0) 

J 2 1 ^ 1 

(j — m + (cosh0) —(cosh0) 

J+ 2 2 


leading to an exact (integral) formula for the fermion 
Casimir energy between the wall and the sphere. 

Asymptotic behaviour. While o-dia) can be used to 
compute numerically the fermion Casimir energy in our 
setup, here we are interested in the asymptotic behaviour 
at large, and, in particular, at short distance, both of 
which can be extracted from the main integral formula 
at the price of some lengthy computations. 

First of all, we can compute the the Casimir interaction 
energy at large separation, i.e., L ^ R. In this regime, 
the dominant contributions are those with j = 1/2 and 
m = ±1/2, and a straightforward computation gives 


Ec 


hcR^ 
ttL^ ’ 


(13) 


leading to an attractive interaction at large distances. 
It seems interesting to notice that the fermion Casimir 
energy, at large distance, falls off as L~^, less rapidly 
than the electromagnetic contribution that decays as 
Compared to the scalar contribution, instead, the 
behaviour in the fermionic case is intermediate between 
that of a scalar with Dirichlet boundary conditions and 
that of a scalar with Neumann boundary conditions, in 
which cases the Casimir energy decays as L~^ and L~'^, 
respectively. 

The more interesting (and computationally more te¬ 
dious) limit is that of the Casimir energy at small sep¬ 
aration, d <C P, beyond the leading (proximity force) 
approximation. In the present case, taking j as the main 
quantum number and using the invariance of the matrix 
Mjm,j'm under the change m i—>■ —m (this implies that 
the next-to-leading order term is of order d smaller than 
the leading term), a lengthy computation returns 


Ecas 


Tn^hcR ( 

2880^2 Y 


1 

3 


20 1 d 


(14) 


It is easy to check that the leading term above coin¬ 
cides with proximity force result, and gives rise to an 
attractive interaction, while the sub-leading term cor¬ 
rects the proximity force result by a positive amount. 
Contrary to what happens at large distance, in this limit 
the fermion and electromagnetic (for perfectly conduct¬ 
ing boundaries) contributions are both attractive and 
scale in the same way. Interestingly, the correction to the 
proximity force result is, for the present setup, positive. 


ddsinh6»e-2«^'^“’^‘'®x 



m ± 1)PY, \ (coshd) 

P/7i=(coshd) 

J + 2 


(/' -I- m)P'^_l (coshd)\ 

-PT^Mcoshd) j’ 
0 2 / 

( 12 ) 


in contrast to the analogous correction for the electro¬ 
magnetic case with perfectly conducting boundaries, for 
which the correction is negative. 

Conclusion. The possibility of manipulating defects in 
condensed matter fermionic systems has triggered new 
curiosity in understanding the analogous of the Casimir 
energy in a fermionic environment and motivated the 
present work. In this paper, we have adapted the mul¬ 
tiple scattering formalism to derive the Casimir interac¬ 
tion energy between a spherical defect and a wall me¬ 
diated by massless fermionic quanta. We have obtained 
an integral representation for the quantum vacuum en¬ 
ergy that is divergence free and valid at all distances. 
From this integral formula, we have extracted the lead¬ 
ing contributions at both large and short distance. The 
behaviour of the Casimir energy at large distances scales 
as L~^ and dominates over the electromagnetic contribu¬ 
tion (for perfectly conducting boundaries), for which the 
energy scales as L~'^. The more interesting result comes 
from the short distance asymptotic behaviour, where the 
leading order contribution is found to coincide with the 
result obtained from the proximity force approximation. 
We have also derived the correction to the proximity force 
approximation, that, in contrast to the electromagnetic 
case, turns out to be positive. As a byproduct, we have 
derived the translation matrices, relating the plane waves 
basis to the spherical wave basis, a result that might be 
useful in other contexts. While the force is attractive 
at both small and large distance for the present choice 
of boundary conditions, it is important to ask how the 
result changes for different boundary conditions (for in¬ 
stance, introducing a phase at one of the boundaries), as 
well as when thermal effects are switched on. Work in 
this direction is in progress. 
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Appendix A: I. Expansions 


A generic solution to the Dirac equation can be expressed as a superposition of the solutions presented in the 
main text (formulae (2) - (5)). In the region between the sphere and the plane, can be represented in two ways. 
In terms of the full set of spherical solutions system (in spherical coordinates centered at O): 

/ OO _ _ 

T. E 

or in terms of the full set of plane waves (in rectangular coordinates centered at O' = Ae^): 

^W(x',t) =H^ 

J — OO j — OO j — OO 

X w) + 4t!2V'ktl2"®(x', w) + + 44.2V'ktlr‘(x', w)) . 

Here x' = x — L, L = Le^. The two representations are related by translation matrices V and W: 




(±).reg^^, 
ki,2 


(x',w)= ^ 




E (42E4S’r®(x, uj) +®(x, w)), 


/(lh),OUt/ N ^2 f f dhy /Txr(=h),ll /(±),OUt/ f \ Trr(lt),21 ,(±),OUt/ / 

(x,w) =H i (x (x , 

/(±),OUt/ \ 7 j2 f f dky /jrT(zL),12 /(±),OUt/ / N . Txr(lt),22 /(±),OUt/ f 

4m.2 (x, w) =H ^ (w^k,w) + w"4:i™4,,2 , 


Using the above expressions, we obtain the following relations 


_„2 r dK r dky 


^ki. 

^k, ,ly 


■ = E 


E 


i=Uf,... 


k_L.im k_Ljim 


f,(±) N 

^lm,2} 


Matching the boundary conditions on the sphere gives 

/k(±) 


1 = -Tit^ I 


°3m,23 


jm I (±) 

U"i.2 / 


while solving the boundary conditions on the plane gives 


(±) 1 42i' 


V^ki,2/ \“kj_,2/ 

These enter the Casimir energy as written in the main text, formula (1): 

h 


ddcas — 


27r 


d^ E Tr In (l - > 


jm,j m 


('IS'r <& 


/ °° dk f 

t;(±)A 2 I -k 


jm,k_i 


kj_ j'm' 


lit 

y(±(.22 

' kj_ j'm' 


where 
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where the following correspondence 

Ti=[T,„,], T2 = [TkJ, = 
is understood. 


foil . 
fo.21 

’ jm.k± 


fol2^ 

fo22. 

' jm,kj 


G^i = 


^kl,/ 


'^k!o'-' 

WP ; , 

kxO m' 


Appendix B: II. Transition Matrices 

The first task is now to use the actual boundary conditions, 

{I + = 0 , 

boundary 

to derive an explicit expressions for the matrices Tjm and Tk^. 

On the exterior of the sphere, using 


(l + iY) 


= 0 , 


r—R 


gives 


(JAkR) T Jj+i{kR)) + T H^/MkR)) = «, 

4^,2Cr {J,+i{kR) ± J,{kR)) + {H^Si^kR) ± H^l\kR)) = 0, 


from which we obtain 


cfs = *-j+^, = 


'Tp(±) _ ( ^3 

I 


0 


0 T' 


(±).2 


(±) 1 ^ Ij^nR) =p iIjj^i{K,R) 

~ Kj{KR)±iKj+i{KR)' 
(±),2 ^ Ij{kR) ± iIj+i{KR) 
Kj{KR)TiKj+i{KRY 

which gives formula (6) in the main text. At the plane, 

^1 — = 0, 


using 


r(=*=) „;,(=*=), o I ^(=*=) „/,(=*=) 


„(±) .,(±).reg. , X _ 


C!iC.r (x',w) + C.2C:2 = L e 


Cli) 

_[_ fciCri+/C2 0'2—fcsCTs 

k 

d(±) 


^ikix-\-ik2y—ik^ zYiujt 


^ikix-\-'ik2y-\-ik3z' 


with 


«-(£ 


""’d£ 


ka = \ k'^ - kj^, 


we get (ui are the Pauli matrices) 


^ ^ ^ ias {ki<7i + fc20'2) 

k k 


^(±) _ _ fo fos (fciui + ^ 20 - 2 ) 

k k 
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from which we can obtain 


or explicitly 


= 


T, 


(±),11 j .(±),12 


.T, 


‘(±).21 


T. 


’(±),22 


= ± 


y/k'^ — k'j_ kia2 — k2(Ti 


ik 


rpW _ rfi22 

'^k_L - -^k_L 


= ± 


y/k^ - kl 
ik 


rpl2 


= T 


i{ki — ik2) 
k 


rp21 


^ i{ki + ik2) 
k 


The above expressions reproduce formula (7) in the main text. 


Appendix C: III. Translation Matrices 


1. Matrix 


The much more tedious task is to find the translation matrices and G^^, which will be explained in the present 
section. In the following we will use 

k = kxBj: + kyBy + k^e^, r = xe^: + yey + ze^, 

with kx = k sin 9 k cos if k iky = k sin 9 k sin ipk,kz = k cos 9 k ■ 

The first step of our procedure consists in defining the following differential operator Vim 


-O ^ + f dx+idy\"'^(m) f dz' 

ij W (^l + m)\\ ik ) ‘ [tkj’ 


Vl-m =1 


121+ 1 {I — m)\ f dx — id, 
47r (I + m)\ \ ik 


P 


(m) 


One has 


Vime 


zk-r 


Ylm{Oki+k)P'^''^i 


Vimjoikr) =i’-ji{kr)Yimi9,ip), 





(kr) = i^h\^\kr)Yim{9,(p), 


= \l 















which allows to express the mode functions given in formula (2) of the main text as 




2j 


j - m 


2j 


V. 


j-i m+i 


±A 


2j + 2 '''3+h^^+i J 


I j — m + 1 


2j + 2 ^•5+2."* 2 


/ j + TO + 1 


P.4 


fSikr), 


Vim .2 — ‘-76 * 


/ j — TO + 1 


_-p 1 , 

2j + 2 


/ j + TO + 1 


P.4 


2j + 2 f+2.™+ 


V 


±1 




/j + TO 


2j 


P 


J-2’^-2 


J - TO 


2j 


P, 




foikr), 


where f* (kr) is defined in formula (3) of the main text. We may now use the following integral representation 


ho{kr) = 


exp(ifcr) 1 f°° f°° ^ika:X+ikyy±iy/k^-kl-klz 


dkr / dk. 


' —00 —00 


ikr 27r 

and express the spinors w), after some calculation, as 


h 

fx \ I fx rxj, IXy 


z^O, 




1 {j-m)l 


47ri V U + m)\ J 

^ikxX+ikyy+iyJ k'^ —k'^—k'^z 


dk„ 


dkye^i'^-i>'^ 


( {j + to)P™ i"" {cos9k)\ 

J 2 

P"^+5(cos0fc)e*‘"'' 

J 2 

±{j-m + 1)P'"T^ (cos Ok) 

2 

V TP™tMcos0fc)P^''y 


h h‘^ — h‘^ — h‘^ 

rx \ / rx fx^ rxy 


Using the definitions of 9k and ifk given at the beginning of the section, one finds 


f {j + (cos6>fc)\ 


P™+^(cos0fc)e*^^ 
J 2 


±{j - TO + l)p”( 1 "" (cos6>fc) 

2 

V (cos 9k)P'^^J 


( 


= {j + m)P {cos9k) 

J 2 


± 


1 \ 

0 

\JkP‘ — kP^ 
k 

ki + ik2 


^—i—/ 


+ p'"\^ {cos9k)P^’‘ 
2 


± 


0 \ 

1 

fci — ik2 


T 
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Recalling that 
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and using the explicit expression for the plane waves (formulae (4) - (5) in the main text), we arrive at 

1 


W. 


(±).ii 


TT 2 {j — Tn) ! 




w. 


(±).21 


2 H^\I U + my. ky/k-^ - k'i 

7r5 / (j — m)! 1 


kxom 2i/2 Y (j + „)! ky/W^ 

fi) 12 fi) 22 

Obtaining Wj)x jm ^kj_ follows from similar steps, leading to 
d±)42 . 


i^ i (cos 0fe)eVfe^-'=i^, 

2 


iTT 2 / (j — m )! 


1 




w, 


U + my.k^W^'' ^,+i 


(±) 2 ? 

Combining all the gives 


21 (±) 

ki.jm I ^(±),21 t^(±),22 I 

V^^kj_,jm 'k±,jm/ 

__^ l ij-my 1 ^(j+m)P™l=(cos6»fe) - m + l)Pj+i ,„_i (cosdfe)^ 

2H^\j {j+my ky/k"^ -k'i^ I P™i^(cos6»fc)e*‘^'“ -iP,+i,„+i (cos6ife)e*‘^'= ) 


as presented in the main text in formula (10). 


2. Matrix 


In 


order to find we may proceed as follows. First of all, we use the integral representation 


sin(fcr) 1 . 

Jo[kr) = ——— = — J difk J dOk sin 

to express the regular spherical solutions as 




[ (j+w)P™i 

J 2 


dS,anS, 


i^(cos6»fe)e<™-5)^A 
2 

P™+^(cos0fe)e*('"+5)‘^'“ 

•? 2 


V 


±{j-m + l)P™“^(cosdfe)e*(™-5)‘^'= 

^p™+5(cos6ife)e<™+5)^'“ 

2 / 


zkr 


/ 


/,(±).reg/ N _ ^reg 1 / (j "l)! f 


dipk / d0fc sin 0k 
Jo 


(j- 


V 


m + l)P;"i^(cos0fe)P(™- 

-p;;tMcos0fc)e<™^ 

±(j + m)P™-^(cos0fe)P(™- 
J 2 

±P'(^tMcos0fc)e<™^ 
•? 2 

































We now introduce the following operators 


w±) (-ir 

jm,l 2 


/7+TO \ 

2 j 


J - W, 


2j 




V 2j + 2 J+i-™+5 


/ j + m+1 

V Y 2j + 2 


/ j j-m+l \ 

J 2j + 2 '^'+5.-™+^ 


^(±) (-1)^ 


/ j + m + l , 

2j + 2 ' 


j+2-m-i 


±1 


j j + m 


2j 


V, 




/ 7 “ 


2j 




that satisfy 


^5w.l •^jm;r®(x,w) 


• V'S’r®(x,w) 

^5w.2 • V'S!r®(x,w) 




x=0 


47r 


j/j?' '-'mm' 5 


=0, 


x=0 


=0 


x=0 


x=0 


= (_l)-J+5* 

47r 


Recalling that 


C;r(=‘'.“)= E 


rn=-j,-j+l,...,j-l,j 




Cr(x',-)= E 
































we can use the above relations for the operators to extract the matrix elements k^ ' 




x=0 


= U - ^ pm ^-^{m-^)vk ^iy/-k^L ^ 

V U + my. 1-2 

- L,a.) 


x=0 


(j - m + 1) (cos0fc) e-<™-5)‘^'=eVfc^-fci^, 


0' +to)! 

y(±).12 _^_ 1 '|-i+| 4 ^p±ii^t-p(±) _^/,(±),reg 




x=0 


=(_l)-™-5^J IZ—(cos6»fe) e-<™+l)‘^<=eVfc"-'=i^, 


(j + to)! •5“2 


x=0 


= (^_iym+^ ^ij iJ_y^pm+^i (cos6»fe) e-<™+^)v’'=eV'="-'=i^. 


(j + to)! J+2 


Combining everything we arrive at 
^(±),ll y(±),n 

k_L ,7m ^k_L ,7m 

r{±),21 ,/±y,22 


^(i)il2 _ I kj_^jm kj_,jm j _ 


jm,ky 


vy-"’’"" K ■ 




(j+to)! 

as reported in formula (9) of the main text. 


(j + TO)P")i^(cos0fc) -P™^^(cos0fe)e-*‘^^ 
i{j - TO + l)P™y = (cos6»fe) iP")t^(cos6>fc)e“*‘^'“ 

J"r 2 2 


Appendix D: IV. Matrix N 


The matrix N can be obtained by combining the translation and transition matrices leading to 




jm,j m 


,=±{-i) 


-m+i , 


U - w)!(j' - to)! 

(j + rn)\{j' + to)! 


/ Ij^kR) =f iljj^i^KR) 
Kj^nR) ± iKjj^i^nR) 


dkj_ k± 


'I'K 


2ky/k^ - kl ^ 


A B 


0 


^2iyyk? —kA L 


0 ^ 

Ij{KR) ± i/j-|-i(kP) 

Kj {kR) =f iKj+i {kR) j 


where 

A ={j + m){j' - TO + 1)P™ 1 = (cos6»fc)P7 /(cosdfe) - P"^i^(cos6lfc)P?)|'Z (cos6»fe), 

J 2 J 2 J 2 d ~r 2 

B =i (^(j + TO)(j' +TO)Pj^i^(cos6»fc)P™_|(cos6>fc) + P™;t"(cos6»fc)Pjr^|(cos6»A;)j , 

c =i ((j - TO + l)(j' - TO + 1)P”^T^ (cos6»fc)P™^ i" (cos 0k) + (cos 0k)p'^^i (cos6>fc)^ , 

D = -{j -m + !)(/ + TO)P”(^^(cos0fe)P?Ji^(cos6>fe) + P™^" (cos6»fc)P'r^Z (cosdfe). 

J“l“ 2 J 2 2 d 2 
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Using the following relation 

(j - m + 1)(/ - m + (cos (cos Ok) + " (cos 9k)P'^^? (cos Ok) 

={j + m){j' + TO)P™i^(cos6lfe)P”(Jf (cos6»fc) + (cos6>fe)P™J^f (cos6>fe), 

J 2 J 2 J 2 J 2 

we can prove that B = C. Then, noticing that 


A B 
C D 


(j+ m)P'" i = (cos6»fe) -P"(t=(cos6>fc)e \ / (/- m + 1)P'“ / (cos 0^) z(j'+ m)P''‘_/(cos 6»fc) 


J-2 \ 3 

,m+i 


^ i(j - m + l)p;(: u (cos 0fe) zP'(:T = (cos 0fc)e-*^'= 

' J"r 2 J' 2 


we can express -A(~^ in terms of 


pr{z) = 


(-1)^ 




jl+rr. 

dz‘+'' 


1 

P™+if (cos6lfc)e*‘^^ 
•? + 2 


-(z^-l)', m>0, 


-iP'^^l (cos6»fe)e*‘^‘^ 
3 2 


as 


7\ 


p-(^) =(-i)™fq:^^r(^), 

_ ( dj (^kR) ^ iIj-\-i (^kR) 


U - fn)\{j' - to)! 
U + m)!(/ + m)! 


KjIkR) ± iKj+i^KR) 

0 


0 ^ 

Ij{KR) ± iIj+i{KR) 
Kj^kR) =f iKj+i{KR) j 


d6lsinh6»e-2''-^“"‘'® 


(j + m)P™ (cosh0) P™'i^(cosh0) \ /(j'— to + 1)P™ ^ (cosh0) {j'+ m)P^_i {cosh0) 


J 2 1 3 2 ^ 

{j — m + (cosh0) —P™ (cosh0) 


i+i 


f+2 

^7+i (coshe) 


J-2 

-P’r\^(cosh6l) / 

•? 2 / 


where we have defined 


1 0\M(1 0 

0 i P 0 -i 


From the above expression, it is straightforward to prove that 


and that 


from which it follows that 


, = M~ ■, , 

jm,j m "> 


M+ = M+ 


hr t°° 

Pcas = -—Re / dKTr In (I - M+) , 

TT Jo 

from which formulae (11) and (12) in the main text can be obtained straightforwardly. 


(Dl) 


Appendix E: V. Small separation asymptotic behaviour 

In order to compute the behaviour of the Casimir energy at short distance, we first expand the logarithm 


OO OO 


He ^ 

^Cas=^^^Re dK J2 Y. ■■■ Y. ^Y,3im---^Y,3om. 


s=0 


and let 


jo = l, ji'-^j + ni, 


£ = —, kR = ui = 
R 


IV^ — 
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which allows us to write 


He ^OO nCO nOO nOO 

Ec. =si i: jn"' y„ ;vCT i ■“ 7.„^ ^ ^ 


Using the following representations 




1 


P™"^fcosh 9) = (y+ "^~ V 

TT ^fc!(j-i-fc)!' 


,(i-i-2fe)e 




p;:t"(cosh 0 )=^^^ 4 ^E 


1 


TT — ' k\ {j — ^ — k)l 


7 27 —1 —2fc • 2 k 2zfm+i)(p 

a^^cos*^ (/psin ipe ^ 


fc=0 

i+^ 






g(i-i-2fe)fl 

g(j+i-2fc)fl r dipcos^i+^-^>^ ^sin^'^ 


pZ\hcosh9) jj+^+^y- y 

J~^2 'TT 


A;=0 


fc! {j + \-k)\ 




d(pcos2j'+i-2'“>sin2'= 


we obtain, after some algebra, 


M~^ 


^ = - —-\/(/ + Ui - to)!(Z + n^+i - m)!(/ + n* + to)!(Z + Ui+i + to)! 

0 \ 


/ Il+ni{uj) - ill+ni + l{uj) 
Kl+mi^) + ^P^Z+rii + llw) 


V 

00 

E 


Il+ni{i^) + Ui+„.+i(w) 
P^z+ni+i(w) — iKi+m +l(‘^)^ 




E i,/i 


A:!(Z + rii + i - fc)! fc'!(/ + n^+i + 5 - fc')! 

dip' COS' 


^{2l-\-ni-\-ni^i—2k—2k')0 


X / d(pcos2'+2’"^-2'=(psin2'=(^e2*'"^ 


'r.n«2i+2"i+i-2fe' gjj^2fe' g2zmv' 


A B 
C D 


where 


A. — ( I -\- Tii -\- -^ — k 


cos (p 
COS(p L 


{l + rii+i -m + 1) + {I + n,+i + m + 1) 


(’ 1 r , 1 ^ 'E 

-i{ip+tp') 

i(ip+ip’) 

\ 2 y cost/?cos<p' 




C =e® cos p cos (^' {I + Hi — m + \) [I + Ui+i — m + 1) e ^ — (^ + + m + 1) (Z + n^+i + m + 1) ^ 


D — ( Z + rii+i + — — k' 


COS p> 
cos ip' L 


(Z + n, - m + 1) + (Z + n, + m + 1) 


We then make the following substitutions 


9 ^ 9 + On, where sinh 9n = 




and, with the aid of a symbolic manipulation program, we expand up to first order in e keeping in mind that Z ^ e~^, 
ni,m ~ E~'^, p,p',9 ~ i/e. For the term involving Bessel functions, we need to make use of the Debye asymptotic 
behaviors 


hiuz) 
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vr]{z) 


y/2TTV (1 + 2:2)1 




/E^iEE A _ EM) 4 _ ^ 
V2z2(1 + 22)H ■■/ 
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where 


and writing 


ri{z) = y/l + + log 


1 + \/l + ’ 


t(z) = 


vT" 


Ul{t) 


t 5t^ 
8 ~ 


^ ^ (U^) 

^l+Ui + liy^') “t" i^l+Ui + liy^') ^l+ni{^^ -j . ■ ^l+ni + l{^') 

Ki+M 


After integration, we finally arrive at 


-^Cas 


7-K^hcR / 

2880 d 2 V 


1 

3 


20 1 d\ 

r) ■ 


(El) 


Appendix F: VI. Proximity force approximation 

The Casimir energy density on a pair of parallel plates separated by a distance d is given by 

II _ Ttt^Hc 

Casl ) 2880#’ 

Hence, the proximity force approximation to the Casimir energy between a sphere and a plate is 


Let 


Then 


pPFA _ 
^Cas — 


[[ dxdysl^ (l- ^R?-x'^-yA 

pR . _, 

= 27 r J dr (^R + d — \/R^ — r^j . 


V = 


R + d- VR"^ - r2 


f(R+d)/d 

Belt =27rd dv{R + d- dv)£l^^ (dv) 

/ OO 

(dv) 


iTT^hcR j‘°° 1 

1440(12 
lir^hcR 
2880(12 ■ 

This coincides with the leading order term we obtain in the previous section. 
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